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Abstract 

We investigate the thermodynamic properties of 5D static and spherically symmetric black holes 
in (i) Einstein-Maxwell-Gauss-Bonnet theory, (ii) Einstein-Maxwell-Gauss-Bonnet theory with neg- 
ative cosmological constant, and in (iii) Einstein- Yang-Mills-Gauss-Bonnet theory. To formulate 
the thermodynamics of these black holes we use the Bekenstein-Hawking entropy relation and, al- 
ternatively, a modified entropy formula which follows from the first law of thermodynamics of black 
holes. The results of both approaches are not equivalent. Using the formalism of geometrothermo- 
dynamics, we introduce in the manifold of equilibrium states a Legendre invariant metric for each 
black hole and for each thermodynamic approach, and show that the thermodynamic curvature 
diverges at those points where the temperature vanishes and the heat capacity diverges. 
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I. INTRODUCTION 

Black holes can be regarded as thermodynamical systems [l, 2] which radiate Hawking 
thermal radiation of temperature proportional to the surface gravity on the horizon of the 
black hole, and Bekenstein-Hawking entropy proportional to its horizon area js, 4|. Indeed, 
these quantities satisfy the first law of black hole thermodynamics 5] . Finding the micro- 
scopic description of this entropy is one of the most challenging questions in theoretical 
physics. The solution to this problem still remains obscure. According to the no-hair the- 
orems of Einstein-Maxwell theory, electro-vacuum black holes are completely described by 
three parameters only: mass M, angular momentum J, and electric charge Q. All these 
parameters define the heat capacity whose sign, in turn, determines the thermodynamical 
stability of the black hole. 

During the last few decades several attempts have been made in order to introduce 
Riemannian geometric concepts in ordinary thermodynamics. First, Weinhold jo] introduced 
on the space of equilibrium states a metric whose components are given as the Hessian of the 
internal thermodynamic energy. Later, Ruppeiner [3] introduced a metric which is defined as 
minus the Hessian of the entropy, and is conformally equivalent to Weinhold's metric, with 
the inverse of the temperature as the conformal factor. The physical meaning of Ruppeiner's 
metric lays in the fluctuation theory of equilibrium thermodynamics. It turns out that the 
second moments of fluctuation are related to the components of Ruppeiner's metric. 

One of the aims of the application of geometry in thermodynamics is to describe phase 
transitions in terms of curvature singularities so that the curvature can be interpreted as a 
measure of thermodynamic interaction. The study of the relation between the phase space 
and the metric structures of the space of equilibrium states led to the result that Weinhold's 
and Ruppeiner's thermodynamic metrics are not invariant under Legendre transformations 
8|. On the other hand, ordinary thermodynamics is invariant with respect to Legendre 
transformations, i. e., the physical properties of a thermodynamic system do not change 
when different thermodynamic potentials are used 9|]. One might then wonder whether 
the use of non Legendre invariant metric structures in ordinary thermodynamics would 
always lead to results that do not depend on the thermodynamic potential. Indeed, several 
examples are known in the literature in which a change of thermodynamic potential leads 
to a modification of the thermodynamic geometry 10|. Some puzzling results arise also in 
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connection with the use of different metrics in the equihbrium space llj, in the sense that 
for the same thermodynamic system the resuhing geometry can be either flat or curved, 
depending on the chosen metric. Recently, in the analysis of Einstein-Gauss-Bonnet (EGB) 
black holes inconsistencies were also found In this work, we will focus on the study of 
EGB black holes and will clarify these inconsistencies. 

Recently, the formalism of geometrothermodynamics (GTD) was developed in order to 
unify in a consistent manner the geometric properties of the phase space and the space of 



equilibrium states [13|. Legendre invariance plays an important role in this formalism. In 
particular, it allows us to derive Legendre invariant metrics for the space of equilibrium 
states. The thermodynamic phase space T is assumed to be coordinatized by the set of 
independent coordinates $,-E",/",a = l,...,n, where $ represents the thermodynamic po- 
tential, and E"- and are the extensive and intensive thermodynamic variables, respectively 
The positive integer n indicates the number of thermodynamic degrees of freedom of the 
black hole configuration. Moreover, the phase space is endowed with the Gibbs one-form 
Q = - dabl^dE^, with 5ab = diag(l, 1), satisfying the condition 6 A (d6)" 7^ 0. Con- 
sider also on T a Riemannian metric G which must be invariant with respect to Legendre 
transformations of the form ($,^",1") (!>, E'^, with $ = 1) - SatE^P, E" = -1°-, 
ja _ Q] impose the Legendre invariance condition in order to incorporate in GTD 
the invariance of ordinary thermodynamics with respect to changes of the thermodynamic 
potential [ol. It turns out that the condition of Legendre invariance is not sufficient to fix 
uniquely the metric G. In fact, Legendre invariance generates the metric 

G = (d$ - ladE'^f + A {UE""!') {xcddE^dl'') , (1) 

where A is an arbitrary real constant, and ^ab and Xab are diagonal constant tensors. Clearly, 
the diagonal components of these tensors can be normalized by rescaling the coordinates E°' 
and J"", and by choosing the multiplicative constant A appropriately. Then, one can express 
^ab and Xab in terms of the usual Euclidean metric, 5ab = diag(l, 1), and the pseudo- 



Euclidean metric, rjab = diag 
final form of these tensors 



— 1,1,...,!). Physical conditions must be invoked to fix the 
15| . Indeed, it turns out that the choice ^ab = ^ab and Xab = ^ab 
i^ab = Sab and Xab = Vab) leads to a metric which describes systems characterized by first 
(second) order phase transitions. Moreover, the choice ^ab = {^ab — Vab) /2 allows us also 
to correctly handle the zero-temperature limit in a geometric way. We see that Legendre 
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invariance leaves free only the signature of Xah^^ The signature, in turn, is fixed by the order 
of the phase transition under consideration l50|. Since in this work we will analyze second 
order phase transitions of EGB black holes 12] , we choose the metric as 



(2) 



The set (T, 6, G) defines a Legendre invariant manifold with a contact Riemannian struc- 
ture. The equilibrium space £^ C T is defined by the map —t- T or, in local coordinates, 
(f : (E^) ($, i?*^, f^), satisfying the condition (p*{Q) = 0, i. e., on S it holds the first 
law of thermodynamics, (i$ = 6abI"'dE'', and the conditions of thermodynamic equilibrium 
ja _ ^o-bg^^gj^b ^ which relate the extensive variables E"" with the intensive ones Then, 
the puUback ip* induces on S, by means oi g = ip*{G), the thermodynamic metric 
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E° 



9$ , 9$ 



VabS' 



be 



dE^dE'^ 



dE^'dE'^ 



(3) 



For the construction of this thermodynamic metric it is only necessary to know explicitly 
the thermodynamic potential in terms of the extensive variables $ = ^(i?"). It is worth 
mentioning that GTD allows us to implement easily different thermodynamic representations 
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22|. 



Since a thermodynamic system is uniquely determined by the fundamental equation 
$ = 9|, the equilibrium space contains the information about the thermodynamic 

geometry of the system. In this connection, two geometric concepts are very important: 
the distance and the curvature of £. The metric g determines the thermodynamic distance 
ds^ = gabdE°'dE^ along the curve which connects two equilibrium states E'^ and E"" + dE"", 



i. e., two points of the equilibrium manifold £ [13|. As has been shown in 16|, the curves 
satisfying the variational principle 5 j ds = Q and the entropy condition determine ther- 
modynamic geodesies, i. e., extremal curves which describe quasi-static thermodynamic 
processes. 

Notice that in general the signature of gab is not fixed and so the equilibrium manifold can 
be either Riemannian or pseudo- Riemannian. This implies that the thermodynamic distance 



ds^ can be either positive, negative or zero. As a consequence, one can show 25|, |26| that 



there exists a casual structure in £ that permits to connect thermodynamically one state 
with another state, but forbids the connection with some other states. The boundary of 
causality is determined by the states satisfying the condition ds^ = 0, i. e., states that can 
be connected by a reversible process. 



Recently, it was pointed out in [23| that the thermodynamic metric gab can be used to 
define in the equihbrium manifold S the probability P{E°') to find a system within the 
interval E"- + dE'^ as 

^(^")-^|giSe.p . (4) 

Then, turns out to be related to the second moments of fiuctuation, in a way that 
resembles the physical interpretation of Ruppeiner's metric in fiuctuation theory. This opens 
the possibility to relate the thermodynamic curvature of S with the correlation length of 
the thermodynamic system described by In this case, the determinant and the principal 
minors of Qab must be positive definite, a property that could be related to the causal 
structure of S mentioned above. A more detailed analysis will be necessary to clarify this 
point. 

Moreover, the curvature of S is interpreted as a measure of the thermodynamic interaction 
of the system with curvature singularities at those points where phase transitions take place 
2^ . This interpretation resembles the role of curvature in gravity and gauge theories, i. e., 
the curvature is a measure of the field interaction with curvature singularities indicating the 
break down of the theory. 

The above interpretation of the thermodynamic curvature in GTD has been proved in 
several ordinary thermodynamic systems and also in the context of black holes. In the case 
of the ideal gas, the curvature vanishes and the thermodynamic geodesies are straight lines. 
This is in accordance with our physical expectations since the ideal gas has no internal 
thermodynamic interaction. In the case of the van der Waals gas, which has a non-trivial 
thermodynamic interaction 9|], the equilibrium space is curved and the curvature singu- 
larities correspond to the points where phase transitions take place. The corresponding 
thermodynamic geodesies are curved and represent quasi-static processes which, under cer- 
tain initial conditions, end at those points where phase transitions greement 
with the well-known relation between geodesic incompleteness and curvature singularities. 
For other ordinary thermodynamic systems, like several generalizations of the ideal gas and 
the Ising model, we obtained similar results. In the context of black holes [l^ , the results 
are compatible with those obtained in ordinary black hole thermodynamics, i. e., black 
holes with thermodynamic interaction possess a curved equilibrium space with singularities 
at those places where phase transitions occur. In the present work, the correctness of this 



interpretation will be shown for black holes in a particular set of higher dimensional gravity 
theories. 

Since a system with thermodynamic interaction can be either stable or unstable, we 
expect that the curvature of £. reproduces this behavior and so it can be different from zero, 
regardless of the stability properties of the system. In fact, we have found so far in GTD only 
two systems with a flat curvature, namely, the ideal gas 13|, |25| and a particular topological 
black hole with fiat horizon in Horava-Lifshitz gravity {28]. Both cases correspond to stable 
systems. In some sense, this can be expected intuitively since in a stable system with 
vanishing curvature there is no thermodynamic interaction that could drive the evolution 
into another state. On the contrary, an unstable system must naturally evolve into a different 
state, a process that demands a non-zero thermodynamic interaction and, consequently, a 
non-fiat curvature. In all the thermodynamic systems we have analyzed in GTD so far, an 
unstable system is always characterized by a non-zero curvature. We will show in this work 
that this is also true in the case of Einstein- Gauss-Bonnet black holes. 

In five dimensions, the most general theory leading to second order field equations for 
the metric is the so-called Einstein- Gauss-Bonnet (EGB) theory, which contains quadratic 
powers of the curvature. The most general action of the EGB theory is obtained by adding 
the Gauss-Bonnet (GB) invariant and the matter Lagrangian Lmatter to the Einstein-Hilbert 
action 

I = kJ d'x^lR + a{R^ - ARf^'R^, + R'^^^'R^^^s) + Lmatter], (5) 

where k, is related to the Newton constant, and a is the Gauss-Bonnet coupling constant. 
GB extensions of General Relativity have been motivated from a string theoretical point of 
view as a version of higher-dimensional gravity, since this sort of modification also appears 
in low energy effective actions of string theory. The pioneering work in this regard belongs to 



Boulware and Deser (29| . They obtained the most general static black hole solutions in EGB 
theory. The GB term has some remarkable features. For instance, in higher dimensions, it 
is the most general quadratic correction which preserves the property that the equations of 



motion involve only second order derivatives of the metric 30 [.However, in 4D, the GB term 



is topological in nature and it does not enter the dynamics 



icslail. 



The GB term is important 



from both physical and geometrical points of view; it naturally arises as the next leadin 



order of the a-expansion of the heterotic superstring theory (a ^ is the string tension) 32| . 



g 



and plays a fundamental role in Chern-Simons gravitational theories 



33|. Several aspects 



of such extended gravity theories have been extensively studied in [34 1. 

In this work, we study the thermodynamics of static spherically symmetric black holes in 
five dimensional Einstein-Maxwell-Gauss-Bonnet (EMGB) with and without cosmological 
constant, and in the Einstein- Yang-Mills-Gauss-Bonnet (EYMGB) theory. We use two dif- 
ferent approaches to formulate black holes thermodynamics. The first one is based upon the 
use of the Bekenstein-Hawking entropy relation, and the second one uses a different formula 
for the entropy which follows from the first law of black hole thermodynamics. We will 
see that the two approaches lead to completely different thermodynamics that have effects 
on the stability properties and the phase transition structure of black holes. The approach 
of GTD is used to show that there exists a Legendre invariant thermodynamic metric for 
the equilibrium space which in all the cases considered here describes correctly the thermo- 
dynamic interaction in terms of the curvature of the equilibrium manifold, and the phase 
transitions and the point with zero temperature in terms of curvature singularities. 

This paper is organized as follows. Section [TTl deals with the thermodynamics and GTD 
of a particular asymptotically de Sitter black hole solution jsSi of EMGB theory. Section UTTl 
is dedicated to the study of an asymptotically anti de Sitter black hole solution of EMGB 
theory with cosmological constant. A black hole spacetime of EYMGB theory is investigated 
in Section IIVI Finally, Section |V] is devoted to conclusions and remarks. Throughout this 
work we use Planck units in which c = G = h = kg = l. 

II. SPHERICALLY SYMMETRIC BLACK HOLE IN EMGB GRAVITY 

In the case of the EGB gravity minimally coupled to the electromagnetic field, the matter 
component of the action ([5]) is given by 



■'matter 



FapF°''^ , FaP — ^P,a — ^q,/3 • (6) 



Spherically symmetric black holes of the EMGB theory have been investigated very in- 
tensively as possible scenarios for the realization of the low energy limit of certain string 



theories. A 



particular solution which contains special black hole spacetime was 



obtained in 35| (see also 36l-l38| ) by using the following 5D static spherically symmetric line 
element 

ds' = -f{r)de + ^^+r^dnl, (7) 



where dfl'^ is the metric of a 3D hypersurface with constant curvature 6k which has the 
exphcit form 



del + sin^ 0iid9l + sin^ Oidel) , 



{k = l) 
{k = -l) 
{k = 0). 



(8) 



del + sinh^ 9i{d9^ + sin^ ^2^^; 
a~^dx'^ 

Here, the coordinate x has the dimension of length while the angular coordinates (6*1, 62) G 
[0,7r] and {63^1^2) G [0,27r]. 

A particular solution is given by the metric function 

.2 



fir) 



4a 



1 ± 



8a{M + 2a\k\) 8aQ^ 



(9) 



^4 3^6 

where the geometric mass M + 2a\k\ is different from that of Einstein gravity for k = ±1. 
This solution is well defined if the expression within the square root is positive definite. 
For the solution (Q of the EMGB theory to describe a black hole it is necessary that the 
condition /(r) = be satisfied. For the special case k = +1, the roots of this equation are 



1 

r± = - 
2 



independently of the value of the coup 




(10) 



ing constant a. It turns out that in some cases these 
radii determine naked singularities 39[. However, the specific case with a < and k = +1 
corresponds to a solution which is asymptotically de Sitter, and represents a black hole with 
an event horizon situated at 




(11) 



provided (;^)^ < |. In this section, we will limit ourselves to the study of this black hole 
spacetime. 

It is interesting to mention that this specific black hole solution is asymptotically de 
Sitter although the cosmological constant does not appear explicitly in the action ([5]). This 



is a particular characteristic of the EGB theory in five dimensions |39|. Moreover, the fact 
that the radius of the event horizon does not depend on the value of the coupling constant 
a leads to interesting thermodynamic consequences. In fact, we will see that if we use 
the Bekenstein-Hawking entropy relationship the thermodynamics of the black hole differs 
completely from the one obtained by using a modified entropy in which the coupling constant 
appears explicitly. 
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A. Analysis with the Bekenstein-Hawking entropy relation 



Since the surface of the event horizon is given by 



/ / sin^ ^ sin = 27rVf^ , (12) 

_ =0 J<i>=o Jip=o 

the Bekenstein-Hawking entropy of the black hole is |^ = ^^^rj^ 40 1. Choosing the 
constants appropriately, the entropy takes the form S = r^, representing the fundamental 
equation that contains all the thermodynamic information. In the mass representation, 
M = M{S,Q), for the black hole solution presented above this fundamental equation can 
be rewritten as 

M = S-^ + -% . (13) 

1 . Thermodynamics 

Using the energy conservation law of the black hole (i. e., dM = TdS + (j)dQ), one obtains 
the temperature and electric potential of the black hole on the event horizon as 

2 3Si-Q^ (14) 

and 

(j) = -^. (15) 

3^1 ^ ^ 

In the positive domain (S'^^'^ > Q'^/3), the temperature increases rapidly as a function of 
the entropy 5* until it reaches its maximum value at S^^^ = 5(5^/3. Then, as the entropy 
increases, the temperature becomes a monotonically decreasing function. This behavior is 
shown in Fig. [H 



According to Davies 4l| , the phase transition structure of the black hole can be derived 
from the behavior of the heat capacity. Strictly speaking, this implies that first we must 
introduce the concept of "heat", say Qheat, for a black hole. This is a complicated question 
that has not been answered so far, in particular, due to the lack of a physically reasonable 
statistical model for black holes j5l| . For this reason, we use here the analogy with ordinary 
thermodynamics as follows. Using the the thermodynamic potential M{S, Q) in which the 
first law of thermodynamics is expressed as dM = TdS + (j)dQ, we define the "heat" through 
the relationship dQheat = TdS so that dM = dQheat + 4>dQ. Then, following the standard 
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FIG. 1: Temperature T, heat capacity Cq and thermodynamic curvature i? of a charged black 
hole in EGB theory as a function of the Bekenstein-Hawking entropy S for Q = 1. The curvature 
singularities are located at the points T = and Cq — t- oo. 

approach of ordinary thermodynamics, we introduce the heat capacity Cq = (^^^)q = 
{w^)o ^^^^ ^^^^ ^^^^ given by 



Cr 



-3S 



3Si 



(16) 



3S-i -5Q\ 

In the physical region with 3^^/3 _ g2 ^ 0, i. e., the region with positive temperature, the 
heat capacity is positive in the interval < 3S^^^ < 5Q^, indicating that the black hole 
is stable in this region. At the maximum value of the temperature which occurs at 
5Q^ = 0, the heat capacity diverges and changes spontaneously its sign from positive to 
negative. This indicates the presence of a second order phase transition which is accompanied 
by a transition into a region of instability (see Fig{T]). 

Usually, when applied to ordinary thermodynamic systems, the above analysis is asso- 
ciated with a particular statistical ensemble, once the internal energy of the system is well 
defined. In the case of black holes, however, there are several possibilities to define the 
internal energy. For this reason, sometimes the potential M{S, Q), or equivalently S{M, Q), 



42], the canonical ensemble with fixed po- 



is associated with the microcanonical ensemble 
tential 43| and, in principle, one could also associate it with the grand canonical ensemble 



because the system is allowed to exchange charge. Here, we will use this last option. 

Using this notation, we can say that the phase transition structure we have found above 
for the EMGB black hole is based upon the analysis of the grand canonical ensemble. On 
the other hand, it is known that the critical points of the heat capacities may depend 
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on the ensemble. In the case of the EMGB black hole we are considering here, from the 
thermodynamic potential M(S', Q) we can define the potentials 



H = M-(f)Q , F = M-TS , G = M -TS -(j)Q , (17) 

that are usually denoted as the enthalpy, the Helmholtz free energy, and the Gibbs free 
energy, respectively. The enthalpy satisfies the first law of thermodynamics, dH = TdS — 
Qd(j), and can be considered as determining the canonical ensemble. Then, if we assume 
again that the "heat" of the black hole is defined by dQheat = TdS, the heat capacity for 
fixed is given by 

C,.f^UrtU-3S. (18) 



dT \dT 

We see that in this ensemble the heat capacity is always negative for S* > 0, indicating that 
the black hole is unstable. In the limit S* — > 0, the thermodynamic description breaks down, 



52| 



because it also implies that if — )■ and M — )■ oo 

Consequently, the phase structure predicted by is completely different from that of Cq. 
In this context, it is interesting to investigate the behavior of the thermodynamic potentials 
at the points of phase transitions. From Eq. fll3p . we can calculate explicitly all the potentials 
and investigate their properties. In Fig. [21 we illustrate the behavior of the potentials near 
the point where the heat capacity Cq diverges. At this point, all the potentials are well 
behaved and only the Helmholtz energy possesses a minimum. This behavior is illustrated 
in the right plot of Fig. |2]that shows the stable minima of the Helmholtz energy for different 
values of the charge. We conclude that the phase transitions corresponding to divergences 
of Cq occur in a region of stability of F. Nevertheless, after the phase transition the black 
hole becomes unstable. 

These are the main features of the thermodynamic behavior of the charged spherically 
symmetric black hole iQ. 



2. Geometrothermodynamics 

For our geometrothermodynamic approach to black hole thermodynamics all what is 
needed is the fundamental equation, M = M{S,Q) as given in Eq. f|T3|) . Then, from the 
general metric (E]) with $ = M and E"" = [S, Q) we obtain the thermodynamic metric of 
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FIG. 2: Left plot: Heat capacities and thermodynamic potentials of the EMGB black hole for a 
fixed value of the charge Q = 1. Right plot: The Helmholtz free energy for different values of the 
charge Q. The dotted curve denotes the location of the minima where the phase transitions occur. 

the equilibrium manifold: 



points where the temperature vanishes. The second singularity corresponds to the roots of 
35^/3 - 5Q^ = 0. According to the expression for the heat capacity Cq given in Eq.f lTGll . 
these are exactly the points where phase transitions take place and the temperature reaches 
its maximum value. It follows that the geometrothermodynamic curvature of the metric 
( IT9|) reproduces correctly the thermodynamic behavior near the points of zero temperature 
as well as near the points of phase transitions (see FigJT]). The above results show that GTD 
correctly describes the thermodynamic properties of the black hole under consideration in 
the grand canonical ensemble in which Cq was considered. 

Notice that the curvature of the metric ( IT9l) cannot reproduce the phase transition struc- 
ture predicted by C^, because it represents a different definition of phase transitions based 
on the use of a different ensemble, namely, the canonical ensemble. Nevertheless, in GTD 




(19) 



The corresponding scalar curvature is given by 

243 ^8/3 
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we can also reproduce the results obtained from the analysis of by using the fact that the 
general thermodynamic metric ([3]) can be applied to any potential in any representation. 
Below we will show this explicitly. 

Consider the canonical ensemble with the thermodynamic potential 

H^M-ct>Q = S'l-i(l--A , (21) 



in which the first law reads dH = TdS — Qdcf). Then, the corresponding temperature can 
be expressed as 

T = \s-'''{l--/) , (22) 

and the heat capacity C^, as mentioned above, is always negative for S > 0, i. e., the 
black hole is unstable. Since the fundamental equation in this ensemble is H = H{S, (p), we 
choose the thermodynamic potential $ = and the independent thermodynamic variables 
E"' = {S, 0) to construct the thermodynamic metric of the equilibrium manifold. Then, 
introducing Eq.( ET]) into the metric (E]), we get 



(23) 



from which we obtain the curvature scalar 



R = ^ . (24) 

54/3 (1_ 3^2)3 

In GTD, we interpret the non-vanishing of the curvature as due to the presence of ther- 
modynamic interaction in the system, regardless of its stability properties. In this case, we 
see that a non-zero curvature describes the thermodynamics of a completely unstable black 
hole. Indeed, from the above expression for the curvature scalar, it follows that there is a 
first singularity for S ^ that corresponds to the limit C,^ — ?■ at which, as mentioned 
above, the thermodynamic description breaks down. The second singularity is located at 
0^ = 4/3 and corresponds to the limit of vanishing temperature. We conclude that the 
curvature singularities of the equilibrium manifold correspond to the locations where the 
thermodynamic description of the black hole in this ensemble breaks down. 

The above results corroborate the well-known fact that the phase transition structure 
depends on the chosen ensemble. In this context, it is interesting to investigate the thermo- 
dynamic properties that follow from the study of the Gibbs function. From Eqs.f ll4l) and 
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(fT5|) . we obtain 



so that 



S 



27T3 



Q 



2(j) 



1 



G = M -TS - (t)Q 



3T^ \~ 4 
4 /^_3 



(25) 



(26) 



27T2 4^ 

Thus, the Gibbs energy is a well-behaved function in the interval 0^ < 4/3 of positive 
temperature. To study this ensemble in GTD, we consider the thermodynamic metric ([3]) 
with $ = G and E"^ = (T, 0), and introduce in the resulting metric the fundamental equation 
( l26l) . Then, we obtain the expression 

2 

dT^ + (4 - 1502)^^02 
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16 



3T2 



l--< 



(27) 



for which the following curvature scalar can be derived 

729 



R 



(28) 



(1-102)^(4-1502)2 

The first curvature singularity appears at 0^ = 4/3, and corresponds to the limit T — 
at which the thermodynamic approach breaks down. The second singularity at 0^ = 4/15, 
which coincides with the point where the second derivative of G vanishes, indicates the 
presence of a second-order phase transition. This result is in accordance with the one 
obtained for the thermodynamic potential M = M{S,Q) in Eq. (l20|) . Indeed, using the 
Eqs. fl25|) . one can show that 35^/'^ — 5Q^ oc 4 — 150^, and this is exactly the term in the 
denominator of (j20l) that vanishes as Cq — ?■ oo. 

The thermodynamic geometry of this black hole was also studied in 35|] using the Rup- 



peiner geometry. It turns out that the Ruppeiner metric is fiat in this case and, consequently, 
cannot reproduce the behavior at the places where phase transitions occur or the tempera- 
ture vanishes. 



B. Geometrothermodynamics with a modified entropy relation 

Usually the entropy of black holes satisfies the so-called area formula, i.e, the black hole 
entropy equals one-quarter of the horizon area. In gravity theories in higher dimensions and 
with higher powered curvature terms, however, the entropy of black holes does not necessarily 
satisfy the area formula and other possibilities can be considered to define entropy. For 
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instance, in 4^ a simple method was suggested to obtain the black hole entropy, by assuming 
that black holes, considered as genuine thermodynamic systems, must obey the first law of 
thermodynamics. That is, we suppose that a black hole solution, parameterized by the mass 
M or, alternatively, by the outer horizon radius th-, and the temperature T, satisfies the first 
law of thermodynamics dM = TdS + fiidQ\ where are the additional charges of the black 
hole and /ij are the corresponding chemical potentials. If the mass and the temperature can 
be calculated by using standard methods, the integration of the first law yields the modified 
entropy formula 



S 



dM 



dr 



H 



So 



(29) 



10 \drH, 

where the additive integration constant 5*0 can be fixed by imposing the condition that the 
entropy goes to zero in the case of an extreme black hole or when the area of the horizon 
vanishes. Notice that in the integration ( 129|) the charges Qi must be considered as constants. 



In 



45 



46| . the modified entropy was computed for an n— dimensional generalization of the 
black hole solution dHl) with the result 



S 



AG 



n-2 



1 + 



2ak{n - 2) 



(30) 



{n — 4:)r'jj 

where Qk denotes the spatial volume element, and a = {n — 2){n — 3)a. In the case n = 5, 
we are considering here, the modified entropy reduces to 



S = r^j + Qkaru 



(31) 



where suitable units were chosen and we set 5*0 = for simplicity. Notice that the contri- 
bution of the correction term vanishes for /c = so that the GB term has no effect on the 
expression for the entropy which reduces in this case to the standard area formula. More- 
over, the modified entropy formula does not contain the additional charges Qi explicitly, but 
only implicitly through the horizon radius r^. So we can assume the validity of the modified 
entropy (13T!) regardless of the nature of the additional charges. 

For the black hole solution ([9]) the modified entropy formula (with k = 1) gives 

3 



S 



2Q_ 



+ , M 



2Q 



3a 




(32) 



In this case the fundamental equation is of the form S = S{M, Q) and cannot be rewritten 
explicitly as M = M{S, Q). This means that for the further analysis we must use the entropy 
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FIG. 3: Modified temperature, heat capacity Cq, and entropy as functions of the specific charge q 
and the coupling parameter a for a fixed value of the mass M = 1. 



representation. This is not a problem for the formalism of GTD which can be applied to any 
arbitrary representation. In fact, for the entropy representation we only need to consider 
the fundamental one-form as 

es = dS-^dM+^dQ, (33) 

so that the thermodynamic potential is $ = S*, the coordinates of the equilibrium manifold 
are -E" = (M, Q), and the equilibrium conditions become 

dS_ _ _dS_ 
T ~ DM ' T ~ dQ ' ^ ' 

From the above expressions one obtains the temperature and electric potential of the black 

hole on the event horizon as 



3 (4a + M + Myi^) ( v^TTg + v/T^) ' 



where q = represents a rescaled specific charge that satisfies the condition < 1. 
Furthermore, to find out the phase transitions structure we must find the points where the 
heat capacity {Cq = TdS/dT\Q = -{dS/dMf /{d^S/dM^)\Q) 



^ 3 VMv/r^(yr+^+v/r^)(4a + M + Myr^)^ 

^ 4 Mil - 3g2) + (4a + M)^Jl-q^- 8a 

diverges. Since the coupling constant a is negative, the temperature function turns out to 
be positive definite only for certain ranges of a and q. In Figj3l we choose a particular range 
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FIG. 4: The heat capacity Cq as a function of the specific charge q for a = —1/4 (left plot) and 
a = —1/10 (right plot). In both cases the mass is M = 1. 

of values of a G [—1/4,0] in which the modified temperature is positive. We also plot the 
modified heat capacity and entropy in the same range of values. Notice that the entropy is 
not positive definite in this interval; however, it is possible to choose the arbitrary constant 
Sq in Eq. fl29|) so that the modified entropy is always positive and the expressions for the 
modified temperature and heat capacity remain unchanged. 

An interesting result of using the modified entropy is that the phase transition structure 
now depends on the value of the coupling constant a. For instance, if we choose it as a = 
— 1/4, the heat capacity is as illustrated in Fig. |l](left plot). In this case, the heat capacity is 
represented by a negative smooth function with no singularities in the interval — 1 < g < 1, 
indicating that the black hole is a completely unstable thermodynamic system with no phase 
transition structure. This behavior changes drastically, if we choose a different value of the 
coupling constant. In fact, FigJH (right plot) illustrates the singular behavior of the heat 
capacity in the case a = —1/10. We can see that at g ~ ±0.82, the black hole undergoes 
a second order phase transition. In the interval —0.82 < q < 0.82, the configuration is 
unstable because the heat capacity is negative. Outside this interval, however, the black 
hole is stable. We conclude that the coupling constant a can induce a second order phase 
transition in an unstable black hole in such a way that the resulting configuration is a stable 
black hole for certain values of the specific charge q. 

Now we investigate this black hole configuration in the context of GTD. As mentioned 
above, the coordinates of the equilibrium manifold are E"" = (M, Q) and the thermody- 
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FIG. 5: Thermodynamic curvature as a function of the rescaled specific charge q. For a = —1/4 
(left plot) the curvature is completely regular, except at g' = ±1. Curvature singularities appear 
for q w 0.82 and a = -1/10 (right plot). 

namic potential $ = is given by means of the fundamental equation ( l32l) . Then, the 
thermodynamic metric for the equilibrium manifold S is given by 



The behavior of the corresponding scalar curvature is shown in Fig|S]for two different values 
of the coupling constant. The plot on the left shows the case a = —1/4 and corresponds 
to the case of an unstable configuration as shown in Fig. H] (left plot). We can see that the 
curvature is represented by a smooth function that is free of singularities in the entire domain 
of the specific charge, except at g = ±1 where the temperature vanishes (see FigJ3]). The plot 
on the right illustrates the behavior for a = —1/10 and shows two curvature singularities at 
q ^ ±0.82 which are the points where the phase transition occurs (see right plot in Figj4]). 
In this case, it is also possible to show that an additional curvature singularity (not plotted) 
exists in the limiting case g — > ±1, indicating the blow up of the approach as T — t- 0. 



9 = 




(38) 
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III. SPHERICALLY SYMMETRIC BLACK HOLE IN EMGB GRAVITY WITH 
COSMOLOGICAL CONSTANT 



In the case of the Einstein-Maxwell- Gauss-Bonnet theory with cosmological constant, the 
matter component of the action is given by 



■'matter 



F^pF'^P - 2A , F^p = Ap^^ - A^^f, , (39) 



where A is the cosmological constant, and Fa/3 represents the electromagnetic Faraday tensor. 
A five dimensional spherically symmetric solution in EMGB gravity with A was obtained 



by Wiltshire 



43], using the metric ansatz ([Tj) with k = +1 and the metric function 



/(r) = IH WlH . (40) 

•'^ ' 4a 4a V r4 Sr^ 3 ^ ^ 



The two parameters M(> 0) and Q are identified as the mass and electric charge of the 
system. The limit of vanishing cosmological constant generates a solution contained in 
Eq.dnj with the minus sign in front of the square root and a redefined mass parameter. 
In this limit, however, the resulting solution does not describe a black hole, but a naked 
singularity. 

In order for the solution ( 140|) to describe a black hole spacetime, it is necessary that the 
expression inside the square root be positive and the function /(r) vanish on the horizon 
radius, i. e., 



^ 8aM 8ag2 4a A ^ . r% rl I 8aM 8aQ^ 4aA ^ 

iH >0, 1 + — -— WlH = 0. (41) 

rji 3r% 3 ' 4a 4aY rj^ 3r% 3 ^ ^ 

Moreover, to guarantee that the mass of the black hole is always positive (see below) we 
must demand that the coupling constant a be positive and the cosmological constant A be 
positive definite. In this section we will limit ourselves to this range of parameters, so that 
the black hole determined by the function (HO!) turns out to be asymptotically anti de Sitter. 



A. Analysis with the Bekenstein-Hawking entropy relation 

The condition /(r//) = implies that 

- 2r^ + 2 (M - 2a) rl-^Q' = 0. (42) 
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Moreover, as we mentioned in Sec. (Ill Ap . with the appropriate choice of units the Bekenstein- 
Hawking entropy of the black hole is given hj S = rjj. Then, the corresponding thermody- 
namic fundamental equation in the mass representation becomes 

M = 2a + S'/' + ^ - ^S'^' . (43) 

Notice that to guarantee the positiveness of the mass, we must assume that a > and 
A < 0. 



1. Thermodynamics 

Using the energy conservation law of the black hole (i. e. dM = TdS + (j)dQ), one obtains 
the temperature and electric potential of the black hole on the event horizon as 

9 ■ 



and 

. 2Q 



3Sl' 



(45) 



Now, for the grand canonical ensemble the heat capacity has the form 



C -2S ^S-s-AS'-Q' \ 

\5Q^-3Sl-ASy ■ ^ ^ 

The expression for the temperature (144|) shows that it is positive only in the range 
2^4/3 _ ^ ^2 ^ q2 Consequently, the heat capacity can take either positive or nega- 
tive values, indicating the possibility of stable and unstable states for the black hole. In 
fact, the expression for the heat capacity exhibits a very rich structure in the range where 
the temperature is positive. In FigJHl a particular range was chosen to show the behavior of 
the temperature and the heat capacity. The plot on the right shows for the particular value 
Q = 1/2 two different phase transitions at S* ~ 0.6 and S ~ 4.9. The first one corresponds 
to a transition from a stable state {Cq > 0) to an unstable state {Cq < 0). The second one 
represents a second order phase transition in which the black hole becomes a stable system 
again. 

Let us now consider the thermodynamic potential 

H = M-(j)Q = 2a + 5^/3 - ^02^2/3 _ :^54/3 _ (47) 
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FIG. 6: Behavior of the temperature and the heat capacity Cq in the range Q G [—1, 1] and 
S € [0, 10] for a fixed value of the cosmological constant A = —1. The right plot shows the details 
of the phase transition structure for the particular charge Q = 1/2. 

for the canonical ensemble in which the first law of thermodynamics reads dH = dQheat—Qd4> 
with dQheat = TdS, as before. In this ensemble, we can define the heat capacity at fixed 
electric potential, i. e. 



whose behavior strongly depends on the value of the cosmological constant. Figure [7] shows 
that for a given value of A, the heat capacity can be either positive or negative with a quite 
complex singularity structure. If we take into account the condition 125*^/^— 90^5*^/^— 4AS' > 
0, which follows from the condition T > 0, we find that stable states {C(f, > 0) are allowed 
for entropies in the range 



If we choose 0^ > 4/3, this condition is always satisfied for S > 0, indicating that stable 
states always exist in this case. Moreover, from the heat capacity (HHl) we see that the roots 
of the equation 125'^/'^ — 90^5^/^ + 4AS' = determine the locations where second order 
phase transitions occur. This behavior is illustrated in Fig. |8l The left plot, where the 
interval has been chosen such that T > is always satisfied, shows a phase transition at 
Sc ~ 0.61. The black hole is stable for S > Sc, and unstable in the interval < S < Sc- 
The right plot shows that for a given value of the entropy, it is possible to find a range of 
values for in which the temperature is positive definite and the heat capacity is positive, 
singular or negative, indicating the existence of stable and unstable black holes. 




12^V3 _ 9^2^1/3 _ 
125V3 - 902^1/3 + ' 



(48) 




(49) 
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FIG. 7: The behavior of the heat capacity as a function of the entropy S for different values of 
the electric potential cp and A = — 1. 
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FIG. 8: The behavior of the heat capacity (dotted curve) and the temperature (solid curve) for a 
fixed value of the electric potential, (/> = 1, (left plot) and for a fixed value of the entropy, S = 2 
(right plot). In both plots we set A = — 1. 

Notice that the singularities of are different from those of Cg; consequently, the 
corresponding phase transition structures do not coincide. In fact, in Fig. [9] we show the 
behavior of the heat capacities and all thermodynamic potentials for this case. First, we 
see that at the points of phase transitions in all the potentials are well behaved and no 
critical points are observed. Moreover, the first divergence of Cq (at S ~ 0.62) is situated 
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FIG. 9: Heat capacities and thermodynamic potentials in terms of the entropy for a = 1, A = — 1, 
and Q = 1/2. For clarity, two different intervals of S are depicted with different scales. 

on a local minimum of F whereas the second singularity (at S* ~ 4.8) is located on a local 
maximum. This situation is also illustrated in Fig. [TO] where the Helmholtz free energy is 
plotted in terms of the entropy for different values of the electric charge. The first phase 
transition occurs on a metastable point of F and describes a configuration in which a stable 
black hole transforms into an unstable one. On the contrary, the second phase transition 
occurs on an unstable value F and corresponds to a black hole that passes from an unstable 
state to a stable state. 



2. Geometrothermodynamics 

In this subsection, we derive the thermodynamic metrics for the equilibrium manifold. In 
the case of the grand canonical ensemble, the fundamental equation is given as M = M{S, Q) 
in Eq. (H3|) . Then, we associate the coordinates E"" = {S,Q) to the equilibrium manifold S 
and $ = M is the thermodynamic potential. The thermodynamic metric ([3]) can then be 
written as 



-^(35^/=^ + AS^ - 5Q^)dS' + dQ' 



(50) 



A straightforward computation results in the following scalar curvature: 



^^27 Sy^NjS, g. A) 

2 (3 54/3 _ g2 _ A ^2)3 (3 54/3 _ 5 g2 + A 52)2 ' 
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FIG. 10: The Helmholtz free energy as a function of the entropy for different values of the electric 
charge. The heat capacity Cq (dotted curve) is also depicted for Q = 0.5 to illustrate that phase 
transitions occur on the minimum and maximum of F for the same value of Q. The left plot covers 
a larger interval of the vertical axis to illustrate the behavior of Cq. The right plot illustrates in 
more detail the behavior of F. 



From the expression for the scalar curvature it is obvious that the singularities are located 



where T — )■ 0, and at the points satisfying the equation 3 S^^'^ — 5 + A S*^ = 0, which 
are the points where Cq — > oo. For instance, for the particular case A = — 1 and Q = 1/2 
the singularities are shown in FigJTT] their locations clearly coincide with the points where 
second order phase transitions occur (see right plot in Figj6]). 

To see if the phase transitions predicted by can also be described in the context of 
GTD, let us consider the thermodynamic metric ([3]) with the the canonical. In this case, the 
thermodynamic potential is $ = and E"" = {S, (p) are the coordinates of the equilibrium 
manifold. Then, using the fundamental equation (147|) . we obtain from ([3]) the metric 



with 



N{S,Q,A) 



42 Q^S^/^A - 34 SQ^A - 5 S^Q^A^ - 18 Q^v^ 
- 7 S^A^ + 36 5"/^A + 15 S^^/'^A^ - 162 + 108 Q^S^/^ . 



(52) 



at the points satisfying the equation 3 S^^^ — — A 5*^ = 0, which coincide with the points 




1<P' - §5^/^) [^S^'^' (l - 1<P' + f 5^/-^) dS' - S'/^'^ , (53) 
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FIG. 11: Curvature singularities in the equilibrium manifold of the black hole (|40p with A = — 1 
and Q = 1/2. The singularities are located at S « 0.6 and S « 4.9. 



and the corresponding curvature scalar 



R 



iV(^,0,A) 



(12^1/3 _ 902^1/3 _ 4A^)3(125V3 - 902^V3 + 4AS) 



2 ' 



(54) 



with 



N{S,(f),A) = (125^/^-9025^/^-4A^)2+A^^/3[2702(4-30^) + 2A5^/^(9(/.2-4A^2/^-4)] . 

(55) 

The curvature singularities situated at the roots of 125*^/^ — 90^5'^/'^ + 4AS' = determine 
the phase transition structure of the black hole, because they coincide with the singularities 
of C(p. The second set of singularities for which 125"^/^ — 90^5*^/^ — 4AS' = corresponds 
to the limit T — )■ and indicates the break down of the thermodjTiamic description of the 
black hole. 



B. Geometrothermodynamics with a modified entropy relation 

The modified entropy relation (131 p . with k = +1, cannot be solved in this case to obtain 
an explicit fundamental equation M = M{Q,S). We must therefore consider the implicit 
fundamental equation determined by the relationships 



a A 
S = rj, + 6arH, M = - + ^ + rj, - -rj, . (56) 

25 



Then, the main thermodynamic variables can then be expressed as 



^ 2(3r|, - - Ar%) 

<^=^> (58) 

^ 3rH(r|, + 2«)^(34-g^-Arf,) 

6aQ2 + 5Q2^2^ + Q^^^ _ 3^6^(1 + 2aA) - Arf^ ' ^ 

For a physically reasonable configuration we demand the positiveness of the temperature; 

this implies that 3rjj — Ar^ > Q^. For a given value of A and Q, this condition determines 

a minimum horizon radius r^*" for which the temperature is positive. Moreover, from the 

expression for the heat capacity (159|) and from the condition of positive temperature, it 

follows that if the condition 

«|A| > ^ (60) 

is satisfied, the heat capacity is positive and, consequently, all possible black hole configu- 
rations are stable. This is an interesting condition that relates two fundamental constants, 
namely, the tension of the string, proportional to a~^, and the cosmological constant A. 

For the range a|A| < 1/2 where unstable states in principle can exist, let us consider 
the parameters A = — 1 and Q = 1. This choice together with the positiveness condition of 
the temperature fix the value of r™" ^ 0.73 (see left plot in Fig JT2l) . Notice that the value 
of r™" does not depend on the value of the coupling constant a. We explore the behavior 
of the heat capacity in FigHSlfor the entire range a G (0, 1/2), according to the condition 
a|A| < 1/2. One can see that the heat capacity is represented by a smooth positive function 
in the entire domain. We conclude that also in this case all the black hole configurations 
are stable. 

We now investigate the geometric properties of the equilibrium manifold. According to 
the implicit fundamental equation, the thermodynamic metric (E]) can be written as 

g = -hijH, Q, «, A)| (5gV^ + 6ar| - 3^(1 + 2aA) - Arf,) dS^ + <9rH (r| + 2af dQ^^, 

(61) 

where 

/i(r^,g,a,A) - — rUrh + 2aY " ^^^^ 
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FIG. 12: Temperature, heat capacity Cg, and thermodynamic curvature for a black hole of the 
EMGB theory with cosmological constant A = — 1, charge Q = 1, and coupling constant a = 1/4 
(for the temperature) and a G (0,0.5) (for the heat capacity). The temperature is positive for 
rn > r^"" ~ 0.73. 

The expression for the scalar is quite cumbersome but it can schematically be represented 
as 

N{rH,A,Q,a) 

(63) 

where N{rH,Q,-^-,<y) is a finite function in the entire domain of definition. From the ex- 
pression for the scalar curvature, the temperature and heat capacity given above, it follows 
that singularities can take place only at those points where T — or C 00. In FigJT2| 
the behavior of the scalar curvature is shown for a particular choice of the parameters. We 
see that a singularity occurs at the point where the temperature vanishes. The singularity 
situated at {rfj + 6a) = corresponds to the limit — )■ which indicates the breakdown 
of the thermodynamic picture of the black hole and, hence of GTD. No other singularities 
exist because the heat capacity is finite in this domain. 



IV. SPHERICALLY SYMMETRIC BLACK HOLE IN EYMGB GRAVITY 

48| and its properties in Einstein- 



In this section we first describe the black hole solution 
Yang-Mills-Gauss-Bonnet (EYMGB) gravity, and then study the geometry of the black hole 
thermodynamics in the subsequent sections. 

The 5D spherically symmetric solution obtained recently by Mazharimousavi and Halisoy 
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48| has the metric 



ds" = -f{r)dt' + J^ + r^dnl , (64) 



where 



Here M is an integration constant to be identified as the mass, and Q is the only non-zero 
gauge charge. The Yang-Mills field 2-form Fl^^F^""^ = QQ'^/r^ represents in this case the 
matter Lagrangian in the general action ([5]). The metric on the unit three sphere dQ"^ is 
given by 

dnl = ^ {de^ + dcj)"^ + dip^-2 cos 9d(j)dilj), (66) 

with 9 G [0,7r], and {(f),ip) G [0,27r]. The event horizon radius rn satisfies the equation 
fifu) = and is given by 

M = rl-2QHnrH . (67) 

This black hole solution of the EYMGB theory is well defined for all r if the Gauss- 
Bonnet coupling parameter a is positive definite. For a < , the spacetime has a curvature 
singularity at the hypersurface r = r^, where is the largest root of /(r) = 0. 



A. Analysis with the Bekenstein-Hawking entropy relation 

In suitable units, the entropy S of the black hole is given by S* = r|^, and A = 27r^r^ 
is the surface area of the event horizon. Then, according to Eq. fl671) . the thermodynamic 
fundamental equation in the M-representation becomes 

M = Sl --Q^nS . (68) 
3 

This is the main relationship from which all the thermodynamic properties of this black hole 
can be derived. 



1. Thermodynamics 

The expressions for the main thermodynamic quantities, namely, the temperature and 
the electric potential are given by 

2 ,92/3 _q2 a 
T = - j^, <l>=--Q\nS. (69) 
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FIG. 13: Temperature and the heat capacity Cq in terms of the Yang-Mihs charge Q. To ilhistrate 
the thermodynamic behavior, the Bekenstein-Hawking entropy was chosen as S" = 10. 



It follows that for the temperature of the black hole to be positive the charge must satisfy the 
condition Q < S^^^. Moreover, for a fixed value of the entropy, the maximum temperature is 
reached at the value Q = 0, indicating that the Yang-Mills charge reduces the temperature 
of the black hole. This behavior is illustrated in FigJT3l 

Now, in the grand canonical ensemble the heat capacity is given by the expression 

Cq = -3^52/3 _3g2 • (70) 

In the region Q < S^^^, where the temperature is positive, the heat capacity diverges at 
those points where Q = S^^^ / \/3, indicating the existence of a second order phase transition. 
In the interval S^^'^/y/S < \Q\ < S^^^, the heat capacity is positive (and T > 0), i. e., the 
black hole configuration is stable in this interval. Furthermore, the heat capacity is negative 
within the interval < |Q| < S^^^/\/3 which corresponds to an unstable thermodynamic 
configuration. Since the heat capacity at Q = is negative, we conclude that the addition 
of a Yang-Mills charge Q to an unstable neutral black hole not only reduces its temperature, 
but also changes its heat capacity until it becomes positive and the system becomes stable, 
if the charge is sufficiently large. The transition from an unstable state to a stable state is 
accompanied by a second order phase transition. This thermodynamic behavior is illustrated 
in FigiH 

Consider now the canonical ensemble that is determined by the thermodynamic potential 

H^M-<PQ = S'/' + l^, (71) 
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FIG. 14: The heat capacity (dotted curve) and the temperature (sohd curve) in terms of the 
entropy. Here cp = 1 and A = — 1. 



from which the dual thermodynamic variables are obtained as 

165^/3 In^ 5-902 



24S\n'S ' " 41n5 
Furthermore, the heat capacity at fixed potential is given by 

n-\^^\ - 351n5(1652/3ln2 5-9(/)2) 
~ df),^ ~ -16^2/3 5 + 2702(2 + In S) 



(72) 



(73) 



and predicts a phase structure different from that of Cq as given in Eq. fl70|) . Indeed, in 
Fig. [14]we illustrate the behavior of and the temperature as functions of the entropy. In 
the interval S G (0, 0.4) with positive temperature, the black hole is unstable. Then, in the 
interval S G (0.4, 1.8) no black holes can exist with positive temperature. In the following 
interval S G (1.8, 3.7) all black holes are stable and have positive temperature. The vertical 
dotted line situated at 5 ~ 3.7 denotes the singularity of that corresponds to a second 
order phase transition during which the stable black hole changes into an unstable state. 

To investigate the position of the phase transitions in the parameter space we consider 
all the thermodynamic potentials M, H, F, and G. Their behavior is depicted in Fig. [151 It 
can be seen that at the point where diverges, all the potentials are well-behaved with no 
extrema. As for the heat capacity Cq, the location of the singularity coincides with the only 
minimum of F. In Fig. [161 the behavior of F and Cq are depicted for different values of 
the charge. It follows that the phase transitions take place in regions where the Helmholtz 
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FIG. 16: The Helmholtz free energy F (solid curves) and the heat capacity Cq (dotted durves) in 
terms of the entropy for different values of the entropy. In each case, the minimum of F coincides 
with the divergence of Cq. 

potential is stable. 
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2. Geometrothermodynamics 



To investigate the geometry of the corresponding equihbrium manifold in the grand 
canonical ensemble, we use the general metric with $ = M and E°- = {S,Q), and 
the fundamental equation M = M{S, Q) as given in Eq. fl68p . Then, we obtain 

9 = ^ {S'^' - Q') [^^^^^^dS' -6\nS dQ^^ . (74) 



The corresponding thermodynamic curvature scalar can be represented as 

R = ^(^-^) (75) 

(52/3 _Q2)3(52/3_3g2)2 ^ ' ^ ^ 

where N{S, Q) is a well-behaved function of its arguments. We see that there are several 
places where true curvature singularities can exist. First, if Q = S^^^ the curvature scalar 
diverges and, as described above, the temperature vanishes. Then, at Q = S^^'^/\/3 there 
exists a singularity whose location coincides with the values at which the heat capacity 
Cq diverges and second order phase transitions occur. Finally, if In S* — )■ the curvature 
scalar diverges. We interpret this additional singularity as related to a second order phase 
transition which is not contained in Cq. In fact, in analogy to the heat capacity Cq defined in 
Eq.dTOD, we can introduce the capacitance Cs = {dQ/d^)s = {d^/dQ)^^ = {d^M/dQ%\ 
Then, from the fundamental equation f l68|) . we obtain Cs = —3/(4 In 5) so that in the 
limit S* — )■ 1 a second order phase transition occurs. This proves the physical origin of the 
additional singularity of the thermodynamic curvature. The behavior of this thermodynamic 
curvature is depicted in FigJTTl 

We now investigate the geometric properties of the equilibrium manifold of this black 
hole in the canonical ensemble. The fundamental equation H = H{S, (p) is given by Eq. fITT]) 
so that we can identify the thermodynamic potential as $ = if and E"- = (S*, 0). Then, 
from the general thermodynamic metric ([3]) we obtain 



1652/3 In^ S - 902 
' = h?5 



-1652/3 In^ S + 2702(2 + In 5) 



(76) 



I 1 ^r. 1.^3 Q , -2,01/3 I 97^2 Ql/3 In ,Q^2nfi,Qln2 Q _ ^ ' 



5452 In' 5 

whose curvature scalar can be represented as 

(-1651n3 5 + 5402^1/3 + 2702^1/3 In5)2(1651n2 S - 9025V3)s 
The first term in the denominator of R coincides with the denominator of C^ and, conse- 
quently, reproduces exactly the phase transition structure of the black hole. The second term 
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FIG. 17: The thermodynamic curvature with 5 = 10 as a function of the Yang-Mills charge Q. 
Singularities are due to the existence of a second order phase transition at Q ~ ±1.24 or to the 
vanishing of the temperature at Q ~ ±2.15. 

in the denominator of R vanishes as T — 0, indicating the break down of the thermodynamic 
and geometrothermo dynamic descriptions. 



B. Geometrothermodynamics with a modified entropy relation 

In this case, the thermodynamic fundamental equation cannot be written explicitly. 
Therefore, we use the implicit equation 5* = S{M, Q) determined by the expressions 

S = rj^ + Garn , and M = rjj - 2Q^ Inrn ■ (78) 

Then, we find the following thermodynamic variables 

(P = 4QlnrH, (80) 
3rH{rj,-Q'){rj, + 2ar 
-rjj + (2a ± 3Q^)rl + 2Q^a ' ^ ^ 

Notice that in this case the condition for a positive definite temperature reads rjj > Q^. 
Moreover, the explicit presence of the coupling constant a in the heat capacity leads to 
the possibility of modifying the phase transition structure of the black hole by changing 
the value of the GB coupling constant. Indeed, the expression for the heat capacity (IHTj) 
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FIG. 18: Temperature, heat capacity, and scalar curvature in terms of the horizon radius rfj of a 
black in the EYMGB theory. Here the modified entropy relation is used with Q = 1 and (5 = 1. 



diverges for 



rl = lQ^ + a±l ^/9Q^ + 20 aQ^ + Aa^ , (82) 



indicating that for a given value of the Yang-Mills charge it is possible to find a range of 
values of a for which second order phase transitions take place. This behavior is illustrated 
in FigM 

For the thermodynamic system determined by the fundamental equation f lTHj) . the Leg- 
endre invariant metric ([3]) is given by 

9 = -^ ^'^l(^,?+2a")5^ { i-'H + (2« + Sqyi + 2Q'a] c/^2+18r|,(r|,+2«)3 In r^dQ^j, 

(83) 

The expression for the scalar curvature can be schematically written as 

R= N{rnQ,a) 

[-r\j + (2a + 3g2)r2^ + 2Q^&f [r^ - Q^irj, + 6a)3(ln mY 

where N{rH,Q,<y) is a function that is finite at those points where the denominator 
vanishes. We see that curvature singularities occur at rjj = Q^, which is the point 
where the temperature vanishes, and at the roots of Eq. (l82|) . which determine the points 
where Cq — oo and second order phase transitions occur. The singularity situated at 
Inrn corresponds to a second order phase transition determined by the capacitance 
Cs = {dQ/d^)s = — l/(41nr//), according to Eq.(178l). Finally, the singularity situated 
at [rf^ + 6a;) = corresponds to the limit 5 — )■ which indicates the breakdown of the 
thermodynamic picture of the black hole and, consequently, of GTD. A particular example 
of the location of these curvature singularities is depicted in FigJTSl 
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V. CONCLUSIONS 



In this work, we analyzed the thermodynamics of static spherically symmetric black holes 
in the five dimensional Einstein- Gauss-Bonnet theory and its generalizations including an 
electromagnetic Maxwell field (EMGB), a cosmological constant (EMGBA), and a Yang- 
Mills field (EYMGB). This kind of black holes was also recently investigated in 49| with 
results which are compatible with the ones obtained in the present analysis. To investigate 
the thermodynamics of these black holes we use two different approaches. The first one 
is based upon the Bekenstein-Hawking entropy relation, according to which the entropy of 
a black hole is proportional to the area of its event horizon. The second approach uses 
as starting point a modified entropy relation that follows from the assumption that black 
holes satisfy the first law of thermodynamics in higher dimensions. The two approaches are 
not equivalent since the corresponding thermodynamic variables exhibit completely different 
behaviors. In particular, we noticed that the thermodynamics of black holes based upon 
the modified entropy formula depends on the value of the coupling constant of the Gauss- 
Bonnet term that appears in the action of the theory. Phase transitions appear that depend 
on the explicit value of the coupling constant, and change the stability properties of the 
black holes. This kind of phase transitions is absent when the Bekenstein-Hawking entropy 
relation is used. 

To study the phase transition structure that follows from the Bekenstein-Hawking entropy 
we used the original Davies' definition, according to which divergences of the heat capacity at 
fixed charge Cq represent second order phase transitions, and an alternative approach based 
on the analysis of the divergences of C^, where is the electric potential dual to Q. For all 
three black holes studied in this work, we showed that the divergences of Cq do not coincide 
with the divergences of and, consequently, the corresponding phase transition structures 
are different. Moreover, we analyzed the behavior of all possible thermodynamic potentials 
in the parameter space S — Q, and found that the divergences of do not correspond 
to any particular point in the parameter space. On the contrary, the divergences of Cq 
are always situated on points where the Helmholtz free energy F possesses an extremum. 
The remaining thermodynamic potentials do not show any special behavior at the singular 
points. To be more specific, in the case of the EMGB black holes we found that Cq predicts 
the existence of stable and unstable states with a phase transition that occurs in a region 
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of stability of F, and during which the black holes undergoes a transition from a stable 
to an unstable state. For this black hole, the heat capacity implies the existence of 
only unstable states with no phase transitions. In the case of the EMGB black hole with 
cosmological constant, Cq predicts two phase transitions. During the first transition the 
black hole goes from a stable state to an unstable one whereas at the second divergence of 
Cq the black hole undergoes a phase transition from an unstable state to a stable state. 
In the parameter space, the first phase transition turns out to be located in a metastable 
region of F, and the second transition is situated in an unstable region of F. The alternative 
heat capacity predicts stable and unstable states with two different divergences. At the 
first divergence, the black hole undergoes a transition from a stable state to an unstable 
state, and at the second divergence it becomes stable again. None of the thermodynamic 
potentials present extrema on the divergences of C^. As for the EYMGB black hole, the 
corresponding heat capacity Cq contains a singularity that describes the transition from a 
stable state to an unstable one. The singularity is located on a stable point of F. Finally, 
the capacity predicts two divergences the first of which corresponds to a transition from 
an unstable to a stable state whereas the second divergence corresponds to a transition to 
an unstable region. 

For all the black holes analyzed in this work, we use the formalism of geometrothermody- 
namics (GTD) to find the geometric properties of the corresponding manifolds of equilibrium 
states. Once the thermodynamic fundamental equation of the black is given, a standard pro- 
cedure of GTD allows us to compute the explicit form of the thermodynamic metric that 
describes the geometric properties of the equilibrium manifold. It turns out that the ther- 
modynamic metrics depend on the entropy relation used to construct the thermodynamics 
of the black holes under consideration. The thermodynamic metrics obtained from the 
Bekenstein-Hawking relation were derived explicitly for the grand canonical ensemble, with 
the thermodynamic potential M = M{S, Q), and for the canonical ensemble, with potential 
H = H{S,(f)). In the first case, the curvature of the equilibrium manifold reproduces the 
thermodynamic behavior of the heat capacity at fixed charge Cq and, in the second case, 
the phase structure of C<^ is reproduced correctly. We also investigated an alternative ther- 
modynamics based upon the modified entropy relation. In this case, we calculated for all 
the black holes the heat capacity Cq, which predicts a non-trivial phase transition structure, 
and showed that GTD represents correctly the phase transitions as curvature singularities. 
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Moreover, curvature singularities also appear at those points where the temperature van- 
ishes, indicating the limit of applicability of black hole thermodynamics and of GTD. 

In the case of the thermodynamics based on the modified entropy relation, in principle, 
it could be possible to perform a similar analysis with different heat capacities and thermo- 
dynamic potentials. The computations, however, cannot be carried out in a similar manner, 
because it is not possible to write down explicitly the intensive variables in terms of the 
extensive ones. This means, for instance, that we cannot write explicitly the enthalpy H in 
terms of S and and so it is not possible to directly compute the thermodynamic metric, 
using the fundamental equation H{S,(f)). Of course, one could express all the derivatives 
with respect to in terms of derivatives with respect to the radius of the horizon Th that, 
in turn, is a function of S and of the implicit function Q{4>). This is not an impossible task, 
but the computations become rather cumbersome; we do not believe that such an analysis 
would shed more light on the properties of GTD that has shown already to correctly describe 
the thermodynamics based on the Bekenstein-Hawking entropy relation. 

We conclude that the formalism of GTD can be used in the EGB theory in five dimensions 
to describe in an invariant manner the thermodynamic properties of black holes in terms of 
geometric concepts, regardless of the entropy relation used to formulate the thermodynamics 
of black holes. In particular, we found that all the black holes we analyzed in the EGB theory 
are characterized by non-flat equilibrium manifolds. This means that all those black holes 
possess an intrinsic non-trivial thermodynamic interaction. Moreover, since we represent 
the thermodynamic interaction by means of the curvature of the equilibrium manifold, the 
points where the heat capacity diverges and, consequently, second order phase transitions 
occur, are represented in GTD by curvature singularities, indicating the limit of applicability 
of the thermodynamic approach to black holes and of GTD. 

We have seen that the phase transition structure of EMG black holes depend upon the 
thermodynamic ensemble, because the corresponding heat capacities have different singular 
behaviors. It is not clear which structure is the right one. Moreover, the investigation 
of thermodynamic potentials does not seem to clarify this question, because no obvious 
relationship was found between the critical points of the potentials and the divergences of 
the heat capacities. Although we found that the divergences of the heat capacity at fixed 
charge always occur at points where the Helmholtz free energy has extrema, it is not clear 
why we should not take into account divergences of heat capacities that are not related to 
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extrema of the potentials. We believe that a classification of phase transitions that would 
take into account the Legendre invariance of ordinary thermodynamics could shed some light 
into this issue. We expect to consider this problem in the near future. 
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